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We study the competition between ferromagnetic double exchange (DE) and nearest-neighbour antiferro- 
magnetic exchange in CMR materials. Towards this end, a single site mean field theory is proposed which 
emphasizes the hopping-mediated nature of the DE contribution. We find that the competition between 
these two exchange interactions leads to ferro- or antiferromagnetic order with incomplete saturation of the 
(sub)lattice magnetization. This conclusion is in contrast to previous results in the literature which find a 
canted spin arrangement under similar circumstances. We attribute this difference to the highly anisotropic 
exchange interactions used elsewhere. The associated experimental implications are discussed. 

PACS numbers: 75.70.Pa, 75.40.Cx, 75.30.Et, 75.10.Lp 



The colossal magnetoresistance (CMR) manganese ox- 
ides have received considerable attention recently |^ . In 
the CMR regime, these materials exhibit ferromagnetism 
which is generally believed to result from a conduction 
electron-mediated double exchange (DE) mechanism 
In addition, there exists evidence which suggests the 
presence of an antiferromagnetic superexchange of com- 
parable scale. Thus, a full understanding of the magnetic 
order in the CMR materials requires a treatment of the 
competition between ferro- and antiferromagnetic inter- 
actions. The aim of the present paper is to address this 
competition and, by doing so, revisit earlier claims in the 
literature ||^ which suggest that spin canting may be the 
most natural means of accommodating these two oppos- 
ing interactions. 

In the literature, previous related calculations [|| have 
been performed for a strongly anisotropic model in which 
the inter- and intra-layer direct exchange constants have 
different signs §. In this paper it is assumed that di- 
rect interactions have everywhere the same (antiferro- 
magnetic) sign and magnitude. This is viewed as more 
appropriate for the Lai-xCaxMnOa perovskite family 
away from the x=0 endpoint, as well as for the layered 
manganates such as La2-2a;Sri_(_2a;Mn207. In view of the 
considerable interest in these layered systems , the 

present calculations address primarily the two dimen- 
sional (2D) lattice; nevertheless, our 2D results are qual- 
itatively representative of the three dimensional case as 
well. In the present situation, there are no ferromag- 
netically locked layers; in this way spin fluctuations are 
enhanced, thereby leading to stronger fluctuations in the 
electronic kinetic energy. Because these fluctuations are 
not fully captured by lowest order (i.e., Hartree-Fock 
like) treatments, in the present paper we introduce a new 
approach to the problem. 

The present theoretical framework is based on a sin- 



gle site mean field theory which emphasizes the hopping- 
mediated nature of the DE-induced ferromagnetic inter- 
action. This is to be contrasted with alternative ap- 
proaches in the literature which implicitly introduce 
a Heisenberg-like exchange interaction to represent these 
DE effects. Our starting point is the standard Hamil- 
tonian derived for the case of infinite Hund's rule 
coupling: 

H=-t,Y: cos ^ {etc, + etc.} ^ . 5,. 
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Here Cj annihilates a fermion on site j, Si represents a 
classical core (localized) spin {S ^ 1), Jaf is the nearest- 
neighbour antiferromagnetic exchange integral, H - ex- 
ternal field, and cos 6ij — S.i ■ Sj/S"^. Throughout this 
paper, we use units in which the bare hopping to, h, ks, 
fiB, and the lattice period are equal to unity. In the first 
term of Eqn. (|l|), we have omitted Berry phase effects 
which are insignificant for a single-site mean field treat- 
ment. 

The random distribution of localized spins leads in 
Eqn. (^ to a highly disordered electronic hopping prob- 
lem. In our mean field approach we focus on a central 
site characterized by hopping b to the surrounding sites. 
This site is embedded in a medium with average hopping 
t, t b ijl^; for clarity these parameters are indicated 
schematically in Fig. 1. The quantities b and t depend 
in a self consistent fashion on the change, Sft, in the free 
energy, associated with the change in hopping matrix el- 
ements t b. This kinetic energy contribution to SQ, 
which can be evaluated following Ref. [O , is given by 
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+9(b~t)-{cp{zo)-^i-Dt)) , 
where the spectral shift function ^(e) is given by 
1 62 



JeffiM) = JaF - M\ 
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ei^(e) — i'{e) J e — 77 
the bound state energy zq < ~Dt is the root of 



z — rj 



(4) 



i/(e) is the density of states, '^{z) = — Tlnjl + cxp[(/i — 
z)/T]\, /X is the chemical potential, and f{z) = {exp[(z — 
/x) /T] + 1}"-'^ . Eqns. (||-^) are valid for a simple lattice in 
any dimensionality D from 1 to 3; the energy integrations 
are performed over the entire conduction band, extending 
from —Dt to Dt. It should be stressed, that it is because 
of the locality of the perturbation (which represents a 
lattice analogue of an s-wave scattering problem) that 
the quantity 50. can be evaluated exactly pl| ]. 

In the ferromagnetic phase at T > 0, the net energy 
cost of a single-spin fluctuation is (in 2D) given by 

50.1 = 5nDE{hi,t,T) + UAF{coii9i2)2 - Hcosai - 

~4Ja_f (cos 6*12)12 + (cos . (5) 

Here, ^12 is the angle formed by the fluctuating spin Si 
with a neighbouring spin S2, and ai is the angle be- 
tween 5*1 and the average magnetization, M (see Fig. 
|l|). We use the notation (...); to represent averaging 
over the Boltzmann probability distribution of spin Si, 
wi oc exp(— (5il//T). It follows that (cos0i2)2 = Mcosai, 
and 



bl EE (cos2(6li2/2))2 = (1 + Mcosai)/2, 



{bl)i =(l + Af2)/2. 



(6) 



The central mean field equation of our formalism is given 
by AI — (cosai)i. 

In the ferro- and antiferromagnetic phases, it is use- 
ful to construct a reference framework with which to 
compare our results. We define J^ffiM) which repre- 
sents an effective M- dependent exchange constant for a 
Heisenberg-like magnet. The appropriate exchange con- 
stant can be deduced by considering small spin fluctua- 
tions (I cos a — M\ ^ 1), which correspond to small fluc- 
tuations in the hopping matrix elements {\t — h\<^t). A 
perturbation expansion of Eqn. (^, then leads to 



5nDE(b,t,T) K -2 



t-b 



ef{e)v{e)de^2{t-b)\E^\ 
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at leading order in T/t, where is the kinetic energy 
of the carriers for t = 1. In the ferromagnetic state, it 
follows from Eqn. (I?]) that 



(8) 



The second term in the above equation represents the 
DE contribution. This term, which is contained in other 
mean field schemes [|[^, increases as M decreases. As 
a consequence, for moderately strong antiferromagnetic 
exchange interactions Je/ / (M) changes sign as M varies 
from to 1. This behaviour has important consequences: 
it leads to a lack of saturation in the low temperature 
magnetization. Typical results for M{T) are plotted in 
Fig. H for these moderately strong exchange interactions 
{\Eq\+2H < 8Jaf < V2\Eo\+ 2H). Here the solid 
line represents the full calculation, while the dashed line 
is obtained using the effective exchange interaction. For 
comparison we plot (dotted line) the magnetization of a 
conventional Heisenberg magnet with the same Tc- 

The lack of saturation seen in Fig. ^ can be under- 
stood as follows. As T decreases, the magnitude of sin- 
gle spin fluctuations also decreases. This leads to an in- 
crease of M, which in turn implies a decrease in |Je//| 
(and thereby a tendency to decrease M). Thus, through 
this self-adjustment of the effective exchange interaction 
(which never becomes large in comparison with T), the 
magnetization fails to reach its proper saturation value, 
Mq = 1 |13|. These self consistent changes in | Je//| lead 
to inadequacies of the effective exchange approximation 
at low T. As may be seen in Fig. ||, the behaviour 
obtained in this approximation differs significantly from 
that found using the full calculation of M{T). This dif- 
ference is due to the fact that when Jeff ~ T is small, 
quadratic terms (in {t — b)/t) dominate the physics |14|. 

The Neel antiferromagnetic state (of the metallic 
phase) can be treated similarly ||l^. It can be shown 
that Neel ordering (which arises for Jaf > '^~^^^\Eo\ in 
zero field) , always exhibits undersaturation of the sublat- 
tice magnetization. This undersaturation (which leads to 
a finite bandwidth) may be viewed as consistent with the 
presumed metallic state. 

Our discussion thus far has not included the canted 
phase first proposed by De Gennes This phase is 
characterized by spin ordering with two equal sub-lattice 
magnetizations m which form an angle 2j between them. 
In the present model, spin canting requires the presence 
of a magnetic field to break the high degeneracy which 
would otherwise occur. This degeneracy is related to the 
fact that the energy of the system depends solely on the 
cosine of the angle which the spin (of sublattice I) 
forms with its nearest neighbours belonging to sublattice 
II. In the context of single site mean field approaches, 
this energy does not change as the spin Si moves along 
any cone around the average direction of spins of sublat- 
tice II. Thus, on average the spin Si will be aligned with 
sublattice II, rather than I. Therefore, in the absence 
of perturbations (caused by next-nearest-neighbour ex- 
change, anisotropy effects, quantum corrections, or small 
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external fields) the canted state is destabilized. Since the 
underlying degeneracy is site- local , its efi^ects will be 
suppressed only when the energy scale of a perturbation 
per individual spin is comparable with that of the ther- 
mal motion of a single spin, that is, with the temperature 
T. 

To characterize the finite field canted state, we use 
the full non-perturbative expression (0) , with appropri- 
ate modifications to Eqns. ^ ||) |17| ]. We obtain two 
coupled self consistent equations, one for the component 
of (Si) parallel to the magnetization of sublattice I, 

— sin 27 (sinai cos/3i)i + cos27 (cosai)i — m , (9) 

and another for the perpendicular component, 

cos 27 (sin ai cos /3i) 1 + sin 27 (cos ai) 1 = . (10) 

where ai and j3i are polar and azimuthal angles of the 
spin 5*1 measured with respect to the co-ordinate system 
that has, as its polar axis, the average direction of the 
spins of sublattice II. We choose /3i = for the spin 
Si lying within the plane containing the two sublattice 
magnetizations . 

The low-T canted state is found to be stable for 
SJaf > \Eq\+ H. The solutions of Eqns. (|-|l^ for 
typical parameters are illustrated in the inset of Fig. |^. 
One can see that, as T — > in the canted phase, the sub- 
lattice magnetization m approaches its proper saturation 
value m = 1. Note that the ferromagnetic (7 = 0) solu- 
tion to the mean field equations is present at > as 
well. In Fig. || (inset), the corresponding magnetization, 
Mpm{T), is represented by the dotted line. However, 
when the canted (7 > 0) solution exists, it corresponds 
to a lower value of the free energy. The canted solu- 
tion branches from the ferromagnetic one at temperature 
T ^ H; at this point the undersaturated ferromagnetic 
state undergoes a second-order spin-fiip transition into 
the low temperature canted state [Q. This lends sup- 
port to the notion that undersaturation is characteristic 
of the generic low temperature regime |l9|] . 

Typical phase diagrams for the DE-superexchange 
magnet in (left) zero and (right) non-zero field are pre- 
sented in Fig. m For to of the order of an eV, our 
choice of parameters corresponds to reasonable values of 
Jaf ~ 300K. In zero field (left), the solid line repre- 
sents the phase boundary between paramagnetic (PM) 
and antiferro- (AFM) or ferromagnetic (FM) metallic 
phases. The ordered phases are undersaturated at low 
T (for slightly smaller Jaf we find a critical value of 
bandfiUing, xi, which divides the saturated, x > xi, 
and undersaturated regimes). At low temperatures and 
small concentrations (in 2D, x < 0.215), the undersat- 
urated AFM state becomes thermodynamically unstable 
{d^/dx < 0), signalling either the onset of a more compli- 
cated spin arrangement or phase separation. The dashed 



line in Fig. ^ corresponds to the boundary of this region 
{^^i/^x = 0). 

The right panel in Fig. |^ shows that in the presence of 
a magnetic field the PM-FM transition is replaced by a 
smooth crossover (dotted line). The spin arrangement of 
the AFM phase becomes non-collinear (flop-phase), and 
has the same symmetry properties as the canted phase 
(CM), which becomes stable at lower T (replacing the 
H = undersaturated FM and AFM phases). The two 
are separated from the PM and FM region by a second- 
order phase transition, which is represented by the solid 
line. At sufficiently small x the latter approaches the 
H ~ Q Neel transition line. The thermodynamic insta- 
bility line (not shown) is only slightly affected by H. 

These calculations have made a number of predictions 
which can be tested experimentally. The layered ma- 
terials La2-2£cSri+2a;Mn207, with x — 0.4, presumably 
lie either within the region where the system should dis- 
play undersaturated ferromagnetic behaviour at low T, 
or on the brink of this region, where thermal fluctuations 
should still be stronger than in a conventional magnet. 
Some measurements of the absolute value of magnetiza- 
tion in X = 0.4 samples support undersaturation pO| , pl| , 
while others do not |2^. We also note that the presence 
of undersaturation in ferro- and antiferromagnetic phases 
may well signal that in reality the system favours more 
complicated (e.g. spin glass- like, cf. Ref. ) spin order- 
ing, that cannot be addressed within a single-site mean 
field theory. 

It is natural to expect that the relative strength of 
the superexchange interaction is even higher in x — 0.3 
compounds, so that this compound may be suitable for 
observing canting under the proper field and temperature 
conditions: T ~ iJ, although experiments have not yet 
been performed in this regime. There have been no obser- 
vations of an ordered canted state (as distinguished from 
possible canting correlations reported in Ref. Q) in the 
layered compounds. This is consistent with our results. 
Finally, we propose that the magnetization dependence 
of the effective exchange constant (available through spin 
wave measurements) should be studied in both 3D and 
2D systems. 
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Ruckenstein. We acknowledge the support of a Univ. 
of Chicago/ Argonne National Lab. collaborative Grant, 
U. S. DOE, Basic Energy Sciences, Contract No. W-31- 
109-ENG-38, and the MRSEC program of the NSF under 
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FIG. 1. Single-spin fluctuation in the ferromagnetic phase. 
The bold arrow represents the average magnetization, and the 
dashed lines correspond to the hopping amplitude b, which 
differs from the background hopping value t (solid lines). 

FIG. 2. Magnetization vs. temperature in the ferromag- 
netic phase at H — 0, X — 0.4, and Jaf ~ 0.06. The solid, 
dashed, and dotted lines correspond to the 2D DE- superex- 
change magnet, effective exchange approximation, and usual 
Heisenberg ferromagnet, respectively. The inset shows the 
behaviour of the sublattice (solid line) and net (dashed line, 
M — m cos 7) magnetizations in the canted state at H — 0.01 
in comparison with the magnetization of the ferromagnetic 
state (dotted line). 



FIG. 3. Phase diagrams of the DE-superexchange magnet 
for Jaf = 0.06 at 77 = (left panel) and 77 = 0.01 (right 
panel), showing the ferro-, antiferro- (flop-phase at 77 > 0), 
paramagnetic, and canted phases (FM, AFM, PM, and CM, 
respectively) . The behaviour of the system is symmetric with 
respect to quarter-filling, x — 0.5. 
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